Application of Amplicon-Based Targeted NGS Technology for Diagnosis of Drug-Resistant Tuberculosis Using FFPE Specimens.

Next-generation sequencing (NGS) enables rapid identification of common and rare drug-resistant genetic variations from tuberculosis (TB) patients' sputum samples and MTB isolates. However, whether this technology is effective for formalin-fixed and paraffin-embedded (FFPE) tissues remains unclear. An amplicon-based targeted NGS sequencing panel was developed to predict susceptibility to 9 antituberculosis drugs, including 3 first-line drugs, by directly detecting FFPE tissues. A total of 178 tissue samples from TB patients who underwent phenotypic drug susceptibility test were retrospectively tested from January 2017 to October 2019 in the Department of Pathology, Beijing Chest Hospital, China. Phenotypic drug susceptibility test results were used as the reference standard. We identified 22 high-quality mutations from 178 FFPE tissue samples, including 15 high+moderate+minimal confidence-level mutations associated with drug resistance (rpoB D435V, S450F/L; KatG S315T; inhA-fabG promoter c-15t; embB G406S, M306V; rpsL K43R, K88R, rrs a1401g, a514c; gyrA D94G/Y/A, A90V), 6 mutations not associated with resistance (rpoB D435Y, H445S, L430P, L452P; embB G406A/D), and one mutation site embB M306I defined as indeterminate. Compared to the phenotypic method, sensitivities (95% CI) for rifampicin, isoniazid, and ethambutol were 96% (79.65-99.90%), 93.55% (78.58-99.21%), and 71.43% (35.24-92.44%), respectively; while for second-line drugs, it varied from 23.53% (9.05-47.77%) for capreomycin to 86.84% (72.20-94.72%) for streptomycin. Specificities for all drugs were satisfactory (>94.51%). Therefore, important pathological FFPE tissue samples, despite partially degraded DNA, can be used as essential specimens for molecular diagnosis of drug resistant TB by amplicon-based targeted NGS technology. IMPORTANCE Amplicon-based targeted NGS technology focuses on a set of gene mutations of known or suspected associations with drug susceptibility in Mycobacterium tuberculosis (MTB). This method offers many benefits, such as low sequencing cost, easy customization, high throughput, shorter testing time and not culture dependent. Formalin-fixed and paraffin-embedded (FFPE) tissues are important pathological specimen in diagnosing tuberculous disease because they are noninfectious and provide excellent preservation of tissue morphology with low storage cost. However, the performance of amplicon-based targeted NGS method on FFPE samples has not been reported yet. Therefore, we evaluated the performance of this method using FFPE samples collected from January 2017 to October 2019 in the Department of Pathology, Beijing Chest Hospital, China. We demonstrate that the amplicon-based targeted NGS method performs excellent on FFPE samples, and it can be applied to pathological diagnosis of drug resistant tuberculosis.

Avaliação do perfil de microRNAs exossomais em portadores de hipercolesterolemia familial / Evaluation of the profile of exosomal microRNAs in patients with familial hypercholesterolemia

A Hipercolesterolemia Familial (HF) é uma doença hereditária do metabolismo lipídico que causa aos portadores alta incidência de aterosclerose prematura. A HF pode ser diagnosticada clínica e geneticamente, entretanto, apenas cerca de 40% podem ter confirmados pelo diagnostico molecular. Assim, outros sistemas de diagnóstico devem ser avaliados. Ultimamente devido a estabilidade em fluidos biológicos, os exossomos circulantes apresentam grande potencial, pois carreiam um número variado de compostos e são considerados veículos de intercomunicação entre os tecidos. Sabe-se que vários RNAs são carreados nos exossomos, incluindo miRNAs, lncRNA e uma variedade de proteínas. Estes componentes podem ser marcadores de diagnóstico para várias doenças inclusive a HF e suas complicações cardiovasculares. Foram utilizadas amostras de exossomos plasmáticos provenientes de 54 pacientes HF sem uso de estatina por, no mínimo, seis semanas, e 38 indivíduos normolipidêmicos para sequenciamento de miRNAs e estudo da proteômica. Os exossomos foram isolados utilizando dois métodos precipitação química e cromatográfica de exclusão de tamanho e caraterizados utilizando: dispersão de luz dinâmica, Western blotting, rastreamento de nanopartículas (NanoSight), imunomarcação e microscopia eletrônica de transmissão. Os miRNAs e proteínas foram extraídos dos exossomos e analisados por sequenciamento de nova geração e espectrometria de massa, respectivamente. Os dados clínicos, biodemográficos e laboratoriais dos pacientes HF e controles indicaram diferenças significativas esperadas entre os grupos, indicando que foram selecionados adequadamente. A caracterização físico-química dos exossomos mostrou resultados com tamanho de ˜90nm e imunorreação positiva para tetraspaninas. O resultado do sequenciamento identificou acima 2000 miRNAs. Os miR-122- 5p e miR-21-5p apresentaram expressão aumentada no grupo HF (log2FC=1,79 e log2FC=1,27, respectivamente), e o miR-122-5p pós normalização em relação ao controle manteve significativo comparados ao controle (p=0,034). A análise comparativa entre exossomos e plasma total mostrou diferença significativa, pois foram identificadas 239 proteínas (p <0,05) diferentes entre exossomos e plasma. Em exossomos, 17 proteínas foram aumentadas e 21 diminuídas em pacientes com HF em comparação com o controle (p <0,05). Destas, seis proteínas foram mais abundantes em HF e sete proteínas foram menos abundantes em exossomos de pacientes com HF em comparação com o controle. A análise de enriquecimento por bioinformática mostrou que a maior parte dessas moléculas (miRNAs e proteínas) foram relacionadas com metabolismo lipídico, dislipidemia, aterosclerose, doença arterial coronariana, adipogênese. Assim, na busca de novos alvos como potenciais biomarcadores de diagnóstico da HF, nossos resultados da análise integrativa entre os miRNAs e as proteínas exossomais abre novas frentes de pesquisa mais bem direcionadas, para a validação desses miRNAs e proteínas exossomais

Familial Hypercholesterolemia (FH) is an inherited disease of lipid metabolism that causes a high incidence of premature atherosclerosis in patients. FH can be diagnosed clinically and genetically, however, only about 40% can be confirmed by molecular diagnosis. Thus, other diagnostic systems should be evaluated. Lately, due to stability in biological fluids, circulating exosomes have great potential, as they carry a varied number of compounds and are considered vehicles of intercommunication between tissues. Several RNAs are known to be carried on exosomes, including miRNAs, lncRNA, and a variety of proteins. These components can be diagnostic markers for several diseases including FH and its cardiovascular complications. Plasma exosome samples from 54 FH patients without statin use for at least six weeks and 38 normolipidemic individuals were used for miRNA sequencing and proteomics studies. Exosomes were isolated using two methods chemical precipitation and size exclusion chromatography and characterized using: dynamic light scattering, Western blotting, nanoparticle tracking (NanoSight), immunostaining and transmission electron microscopy. MiRNAs and proteins were extracted from exosomes and analyzed by next-generation sequencing and mass spectrometry, respectively. Clinical, biodemographic and laboratory data of FH patients and controls indicated significant expected differences between the groups, indicating that they were appropriately selected. The physicochemical characterization of exosomes showed results with a size of ˜90nm and positive immunoreaction for tetraspanins. The sequencing result identified above 2000 miRNAs. miR-122-5p and miR-21-5p showed increased expression in the FH group (log2FC=1.79 and log2FC=1.27, respectively), and miR122-5p after normalization in relation to the control remained significant compared to the control (p=0.034). The comparative analysis between exosomes and total plasma showed a significant difference, as 239 different proteins (p < 0.05) were identified between exosome and plasma. In exosomes, 17 proteins were increased and 21 decreased in FH patients compared to control (p < 0.05). Of these, six proteins were more abundant in FH and seven proteins were less abundant in exosomes from patients with FH compared to the control. Bioinformatics enrichment analysis showed that most of these molecules (miRNAs and proteins) were related to lipid metabolism, dyslipidemia, atherosclerosis, coronary artery disease, adipogenesis. Thus, in the search for new targets as potential diagnostic biomarkers of FH, our results of the integrative analysis between miRNAs and exosomal proteins opens new and better-directed research fronts for the validation of these miRNAs and exosomal proteins

Genomic prediction using low-coverage portable Nanopore sequencing.

Most traits in livestock, crops and humans are polygenic, that is, a large number of loci contribute to genetic variation. Effects at these loci lie along a continuum ranging from common low-effect to rare high-effect variants that cumulatively contribute to the overall phenotype. Statistical methods to calculate the effect of these loci have been developed and can be used to predict phenotypes in new individuals. In agriculture, these methods are used to select superior individuals using genomic breeding values; in humans these methods are used to quantitatively measure an individual's disease risk, termed polygenic risk scores. Both fields typically use SNP array genotypes for the analysis. Recently, genotyping-by-sequencing has become popular, due to lower cost and greater genome coverage (including structural variants). Oxford Nanopore Technologies' (ONT) portable sequencers have the potential to combine the benefits genotyping-by-sequencing with portability and decreased turn-around time. This introduces the potential for in-house clinical genetic disease risk screening in humans or calculating genomic breeding values on-farm in agriculture. Here we demonstrate the potential of the later by calculating genomic breeding values for four traits in cattle using low-coverage ONT sequence data and comparing these breeding values to breeding values calculated from SNP arrays. At sequencing coverages between 2X and 4X the correlation between ONT breeding values and SNP array-based breeding values was > 0.92 when imputation was used and > 0.88 when no imputation was used. With an average sequencing coverage of 0.5x the correlation between the two methods was between 0.85 and 0.92 using imputation, depending on the trait. This suggests that ONT sequencing has potential for in clinic or on-farm genomic prediction, however, further work to validate these findings in a larger population still remains.

The ENmix DNA methylation analysis pipeline for Illumina BeadChip and comparisons with seven other preprocessing pipelines.

BACKGROUND: Illumina DNA methylation arrays are high-throughput platforms for cost-effective genome-wide profiling of individual CpGs. Experimental and technical factors introduce appreciable measurement variation, some of which can be mitigated by careful "preprocessing" of raw data. METHODS: Here we describe the ENmix preprocessing pipeline and compare it to a set of seven published alternative pipelines (ChAMP, Illumina, SWAN, Funnorm, Noob, wateRmelon, and RnBeads). We use two large sets of duplicate sample measurements with 450 K and EPIC arrays, along with mixtures of isogenic methylated and unmethylated cell line DNA to compare raw data and that preprocessed via different pipelines. RESULTS: Our evaluations show that the ENmix pipeline performs the best with significantly higher correlation and lower absolute difference between duplicate pairs, higher intraclass correlation coefficients (ICC) and smaller deviations from expected methylation level in mixture experiments. In addition to the pipeline function, ENmix software provides an integrated set of functions for reading in raw data files from mouse and human arrays, quality control, data preprocessing, visualization, detection of differentially methylated regions (DMRs), estimation of cell type proportions, and calculation of methylation age clocks. ENmix is computationally efficient, flexible and allows parallel computing. To facilitate further evaluations, we make all datasets and evaluation code publicly available. CONCLUSION: Careful selection of robust data preprocessing methods is critical for DNA methylation array studies. ENmix outperformed other pipelines in our evaluations to minimize experimental variation and to improve data quality and study power.

A scalable high-throughput targeted next-generation sequencing assay for comprehensive genomic profiling of solid tumors.

Timely and accurate identification of molecular alterations in solid tumors is essential for proper management of patients with advanced cancers. This has created a need for rapid, scalable comprehensive genomic profiling (CGP) systems that detect an increasing number of therapeutically-relevant variant types and molecular signatures. In this study, we assessed the analytical performance of the TruSight Oncology 500 High-Throughput assay for detection of somatic alterations from formalin-fixed paraffin-embedded tissue specimens. In parallel, we developed supporting software and automated sample preparation systems designed to process up to 70 clinical samples in a single NovaSeq 6000TM sequencing run with a turnaround time of <7 days from specimen receipt to report. The results demonstrate that the scalable assay accurately and reproducibly detects small variants, copy number alterations, microsatellite instability (MSI) and tumor mutational burden (TMB) from 40ng DNA, and multiple gene fusions, including known and unknown partners and splice variants from 20ng RNA. 717 tumor samples and reference materials with previously known alterations in 96 cancer-related genes were sequenced to evaluate assay performance. All variant classes were reliably detected at consistent and reportable variant allele percentages with >99% overall accuracy and precision. Our results demonstrate that the high-throughput CGP assay is a reliable method for accurate detection of molecular alterations in support of precision therapeutics in oncology. The supporting systems and scalable workflow allow for efficient interpretation and prompt reporting of hundreds of patient cancer genomes per week with excellent analytical performance.

High-Throughput Selection and Characterisation of Aptamers on Optical Next-Generation Sequencers.

Aptamers feature a number of advantages, compared to antibodies. However, their application has been limited so far, mainly because of the complex selection process. 'High-throughput sequencing fluorescent ligand interaction profiling' (HiTS-FLIP) significantly increases the selection efficiency and is consequently a very powerful and versatile technology for the selection of high-performance aptamers. It is the first experiment to allow the direct and quantitative measurement of the affinity and specificity of millions of aptamers simultaneously by harnessing the potential of optical next-generation sequencing platforms to perform fluorescence-based binding assays on the clusters displayed on the flow cells and determining their sequence and position in regular high-throughput sequencing. Many variants of the experiment have been developed that allow automation and in situ conversion of DNA clusters into base-modified DNA, RNA, peptides, and even proteins. In addition, the information from mutational assays, performed with HiTS-FLIP, provides deep insights into the relationship between the sequence, structure, and function of aptamers. This enables a detailed understanding of the sequence-specific rules that determine affinity, and thus, supports the evolution of aptamers. Current variants of the HiTS-FLIP experiment and its application in the field of aptamer selection, characterisation, and optimisation are presented in this review.

High concordance in preimplantation genetic testing for aneuploidy between automatic identification via Ion S5 and manual identification via Miseq.

The Ion S5 (Thermo Fisher Scientific) and Miseq (Illumina) NGS systems are both widely used in the clinical laboratories conducting PGT-A. Each system employs discrepant library preparation steps, sequencing principles, and data processing algorithms. The automatic interpretation via Ion Reporter software (Thermo Fisher Scientific) and the manual interpretation via BlueFuse Multi software (Illumina) for chromosomal copy number variation (CNV) represent very different reporting approaches. Thus, it is intriguing to compare their ability of ploidy detection as PGT-A/NGS system. In the present study, four aneuploid cell lines were individually mixed with a diploid cell line at different aneuploid ratios of 0% (0:5), 10% (1:9), 20% (1:4), 40% (2:3), 50% (3:3), 60% (3:2), 80% (4:1) and 100% (5:0) to assess the sensitivity and specificity for whole chromosomal and segmental aneuploidy detection. The clinical biopsies of 107 blastocysts from 46 IVF/PGT-A cycles recruited between December 2019 and February 2020 were used to calculate the concordance. Initially, the pre-amplified products were divided into two aliquots for different library preparation procedures of each system. Applying the same calling criteria, automatic identification was achieved through the Ion Reporter, while well-trained technicians manually identified each sample through the BlueFuse Multi. The results displayed that both systems reliably distinguished chromosomal CNV of the mixtures with at least 10% aneuploidy from karyotypically normal samples ([Ion S5] whole-chromosomal duplication: 2.14 vs. 2.05, p value = 0.009, segmental deletion: 1.88 vs. 2.05, p value = 0.003; [Miseq] whole-chromosomal duplication: 2.12 vs. 2.03, p value = 0.047, segmental deletion: 1.82 vs. 2.03, p value = 0.002). The sensitivity and specificity were comparable between the Ion S5 and Miseq ([sensitivity] 93% vs. 90%, p = 0.78; [specificity] 100% vs. 100%, p value = 1.0). In the 107 clinical biopsies, three displayed chaotic patterns (2.8%), which could not be interpreted for the ploidy. The ploidy concordance was 99.04% (103/104) per embryo and 99.47% (2265/2277) per chromosome pair. Since their ability of detection were proven to be similar, the automatic identification in Ion S5 system presents comparatively faster and more standardized performance.

Analyzing Modern Biomolecules: The Revolution of Nucleic-Acid Sequencing - Review.

Recent developments have revolutionized the study of biomolecules. Among them are molecular markers, amplification and sequencing of nucleic acids. The latter is classified into three generations. The first allows to sequence small DNA fragments. The second one increases throughput, reducing turnaround and pricing, and is therefore more convenient to sequence full genomes and transcriptomes. The third generation is currently pushing technology to its limits, being able to sequence single molecules, without previous amplification, which was previously impossible. Besides, this represents a new revolution, allowing researchers to directly sequence RNA without previous retrotranscription. These technologies are having a significant impact on different areas, such as medicine, agronomy, ecology and biotechnology. Additionally, the study of biomolecules is revealing interesting evolutionary information. That includes deciphering what makes us human, including phenomena like non-coding RNA expansion. All this is redefining the concept of gene and transcript. Basic analyses and applications are now facilitated with new genome editing tools, such as CRISPR. All these developments, in general, and nucleic-acid sequencing, in particular, are opening a new exciting era of biomolecule analyses and applications, including personalized medicine, and diagnosis and prevention of diseases for humans and other animals.

Real-time next-generation sequencing on shell-vial culture to contribute to diagnosis of lymphatic tuberculosis: a case report.

Lymph node tuberculosis is a of limited clinical suspicion form of Mycobacterium tuberculosis infection. After 15 days incubation in a cellular culture and directly from the supernatant, 11 minutes of Oxford Nanopore MinION sequencing provided a preliminary result of an antibiotic-susceptible M. tuberculosis Indo-Oceanic lineage strain. Oxford Nanopore MinION sequencing is a promising tool for optimising the laboratory diagnosis of lymph node tuberculosis.

Implementation of whole genome sequencing for tuberculosis diagnostics in a low-middle income, high MDR-TB burden country.

Whole genome sequencing (WGS) is revolutionary for diagnostics of TB and its mutations associated with drug-resistances, but its uptake in low- and middle-income countries is hindered by concerns of implementation feasibility. Here, we provide a proof of concept for its successful implementation in such a setting. WGS was implemented in the Kyrgyz Republic. We estimated needs of up to 55 TB-WGS per week and chose the MiSeq platform (Illumina, USA) because of its capacity of up to 60 TB-WGS per week. The project's timeline was completed in 93-weeks. Costs of large equipment and accompanying costs were 222,065 USD and 8462 USD, respectively. The first 174 WGS costed 277 USD per sequence, but this was skewed by training inefficiencies. Based on real prices and presuming optimal utilization of WGS capacities, WGS costs could drop to 167 and 141 USD per WGS using MiSeq Reagent Kits v2 (500-cycles) and v3 (600-cycles), respectively. Five trainings were required to prepare the staff for autonomous WGS which cost 48,250 USD. External assessment confirmed excellent performance of WGS by the Kyrgyz laboratory in an interlaboratory comparison of 30 M. tuberculosis genomes showing complete agreeance of results.

A New Method for Sequencing the Mitochondrial Genome by Using Long Read Technology.

We describe a protocol to prepare a multiplexed mtDNA library from a blood sample for performing a long read sequencing of the mitochondrial genome. All steps are carefully described to get a high enrichment of mtDNA relative to total DNA extracted from the blood sample. The obtained mutiplexed library allows the production of long sequence mtDNA reads up to 16.5 kbp with a quality enabling variant-calling by using a portable sequencer (MinION, Oxford Nanopore Technologies).

iSeq 100 for metagenomic pathogen screening in ticks.

BACKGROUND: Ticks are blood-sucking ectoparasites that play a pivotal role in the transmission of various pathogens to humans and animals. In Korea, Haemaphysalis longicornis is the predominant tick species and is recognized as the vector of pathogens causing various diseases such as babesiosis, borreliosis, rickettsiosis, and severe fever with thrombocytopenia syndrome. METHODS: In this study, the targeted high-throughput sequencing of the 16S rRNA V4 region was performed using the state-of-the-art sequencing instrument, iSeq 100, to screen bacterial pathogens in H. longicornis, and the findings were compared with those using conventional PCR with specific primers. Microbiome analyses were performed with EzBioCloud, a commercially available ChunLab bioinformatics cloud platform. ANOVA-Like Differential Expression tool (ALDEx2) was used for differential abundance analysis. RESULTS: Rickettsia spp. were detected in 16 out of 37 samples using iSeq 100, and this was confirmed using a PCR assay. In the phylogenetic analysis using gltA and ompA sequences of the detected Rickettsia, the highest sequence similarity was found with 'Candidatus Rickettsia jingxinensis' isolate Xian-Hl-79, 'Ca. R. jingxinensis' isolate F18, and 'Ca. R. longicornii' isolate ROK-HL727. In the microbiome study, Coxiella AB001519, a known tick symbiont, was detected in all 37 tick samples. Actinomycetospora chiangmaiensis was more abundant in Rickettsia-positive samples than in Rickettsia-negative samples. CONCLUSIONS: In this study, iSeq 100 was used to investigate the microbiome of H. longicornis, and the potentially pathogenic Rickettsia strain was detected in 16 out of 37 ticks. We believe that this approach will aid in large-scale pathogen screening of arthropods to be used in vector-borne disease control programs.

The Illumina Sequencing Protocol and the NovaSeq 6000 System.

The NovaSeq 6000 is a sequencing platform from Illumina that enables the sequencing of short reads with an output up to 6 Tb. The NovaSeq 6000 uses the typical Illumina sequencing workflow based on library preparation, cluster generation by in situ amplification, and sequencing by synthesis. Flexibility is one of the major features of the NovaSeq 6000. Several types of sequencing kits coupled with dual flow cell mode enable high scalability of sequencing outputs to match a wide range of applications from complete genome sequencing to metagenomics analysis. In this chapter, after explaining how to assemble a normalized pool of libraries for sequencing, we will describe the experimental steps required to run the pools on the NovaSeq 6000 platform.

Integrative reconstruction of cancer genome karyotypes using InfoGenomeR.

Annotation of structural variations (SVs) and base-level karyotyping in cancer cells remains challenging. Here, we present Integrative Framework for Genome Reconstruction (InfoGenomeR)-a graph-based framework that can reconstruct individual SVs into karyotypes based on whole-genome sequencing data, by integrating SVs, total copy number alterations, allele-specific copy numbers, and haplotype information. Using whole-genome sequencing data sets of patients with breast cancer, glioblastoma multiforme, and ovarian cancer, we demonstrate the analytical potential of InfoGenomeR. We identify recurrent derivative chromosomes derived from chromosomes 11 and 17 in breast cancer samples, with homogeneously staining regions for CCND1 and ERBB2, and double minutes and breakage-fusion-bridge cycles in glioblastoma multiforme and ovarian cancer samples, respectively. Moreover, we show that InfoGenomeR can discriminate private and shared SVs between primary and metastatic cancer sites that could contribute to tumour evolution. These findings indicate that InfoGenomeR can guide targeted therapies by unravelling cancer-specific SVs on a genome-wide scale.

Evaluation of Promega PowerSeq™ Auto/Y systems prototype on an admixed sample of Rio de Janeiro, Brazil: Population data, sensitivity, stutter and mixture studies.

Forensic DNA typing typically relies on the length-based (LB) separation of PCR products containing short tandem repeat loci (STRs). Massively parallel sequencing (MPS) elucidates an additional level of STR motif and flanking region variation. Also, MPS enables simultaneous analysis of different marker-types - autosomal STRs, SNPs for lineage and identification purposes, reducing both the amount of sample used and the turn-around-time of analysis. Therefore, MPS methodologies are being considered as an additional tool in forensic genetic casework. The PowerSeq™ Auto/Y System (Promega Corp), a multiplex forensic kit for MPS, enables analysis of the 22 autosomal STR markers (plus Amelogenin) from the PowerPlex® Fusion 6C kit and 23 Y-STR markers from the PowerPlex® Y23 kit. Population data were generated from 140 individuals from an admixed sample from Rio de Janeiro, Brazil. All samples were processed according to the manufacturers' recommended protocols. Raw data (FastQ) were generated for each indexed sample and analyzed using STRait Razor v2s and PowerSeqv2.config file. The subsequent population data showed the largest increase in expected heterozygosity (23%), from LB to sequence-based (SB) analyses at the D5S818 locus. Unreported allele was found at the D21S11 locus. The random match probability across all loci decreased from 5.9 × 10-28 to 7.6 × 10-33. Sensitivity studies using 1, 0.25, 0.062 and 0.016 ng of DNA input were analyzed in triplicate. Full Y-STR profiles were detected in all samples, and no autosomal allele drop-out was observed with 62 pg of input DNA. For mixture studies, 1 ng of genomic DNA from a male and female sample at 1:1, 1:4, 1:9, 1:19 and 1:49 proportions were analyzed in triplicate. Clearly resolvable alleles (i.e., no stacking or shared alleles) were obtained at a 1:19 male to female contributor ratio. The minus one stutter (-1) increased with the longest uninterrupted stretch (LUS) allele size reads and according to simple or compound/complex repeats. The haplotype-specific stutter rates add more information for mixed samples interpretation. These data support the use of the PowerSeqTM Auto/Y systems prototype kit (22 autosomal STR loci, 23 Y-STR loci and Amelogenin) for forensic genetics applications.

Next-generation sequencing technologies: An overview.

Since the days of Sanger sequencing, next-generation sequencing technologies have significantly evolved to provide increased data output, efficiencies, and applications. These next generations of technologies can be categorized based on read length. This review provides an overview of these technologies as two paradigms: short-read, or "second-generation," technologies, and long-read, or "third-generation," technologies. Herein, short-read sequencing approaches are represented by the most prevalent technologies, Illumina and Ion Torrent, and long-read sequencing approaches are represented by Pacific Biosciences and Oxford Nanopore technologies. All technologies are reviewed along with reported advantages and disadvantages. Until recently, short-read sequencing was thought to provide high accuracy limited by read-length, while long-read technologies afforded much longer read-lengths at the expense of accuracy. Emerging developments for third-generation technologies hold promise for the next wave of sequencing evolution, with the co-existence of longer read lengths and high accuracy.

Developmental validation of Oxford Nanopore Technology MinION sequence data and the NGSpeciesID bioinformatic pipeline for forensic genetic species identification.

Species identification of non-human biological evidence through DNA nucleotide sequencing is routinely used for forensic genetic analysis to support law enforcement. The gold standard for forensic genetics is conventional Sanger sequencing; however, this is gradually being replaced by high-throughput sequencing (HTS) approaches which can generate millions of individual reads in a single experiment. HTS sequencing, which now dominates molecular biology research, has already been demonstrated for use in a number of forensic genetic analysis applications, including species identification. However, the generation of HTS data to date requires expensive equipment and is cost-effective only when large numbers of samples are analysed simultaneously. The Oxford Nanopore Technologies (ONT) MinION™ is an affordable and small footprint DNA sequencing device with the potential to quickly deliver reliable and cost effective data. However, there has been no formal validation of forensic species identification using high-throughput (deep read) sequence data from the MinION making it currently impractical for many wildlife forensic end-users. Here, we present a MinION deep read sequence data validation study for species identification. First, we tested whether the clustering-based bioinformatics pipeline NGSpeciesID can be used to generate an accurate consensus sequence for species identification. Second, we systematically evaluated the read variation distribution around the generated consensus sequences to understand what confidence we have in the accuracy of the resulting consensus sequence and to determine how to interpret individual sample results. Finally, we investigated the impact of differences between the MinION consensus and Sanger control sequences on correct species identification to understand the ability and accuracy of the MinION consensus sequence to differentiate the true species from the next most similar species. This validation study establishes that ONT MinION sequence data used in conjunction with the NGSpeciesID pipeline can produce consensus DNA sequences of sufficient accuracy for forensic genetic species identification.

Parallel sequencing of 87 STR and 294 SNP markers using the prototype of the SifaMPS panel on the MiSeq FGx™ system.

Massively parallel sequencing (MPS), or next generation sequencing (NGS), is a promising methodology for the detection of short tandem repeats (STRs) and single nucleotide polymorphisms (SNPs) in forensic genetics. Here, the prototype SifaMPS Panel is designed to simultaneously target 87 STRs and 294 SNPs with forensic interest in a single multiplex in conjunction with the TruSeq™ Custom Amplicon workflow and MiSeq FGx™ System. Two in-house python scripts are adopted for the fastq-to-genotype interpretation of MPS data concerning STR and SNP, respectively. In the present study, by sequencing 50 Chinese Hans and many other DNA samples involved in validation studies, system parameters including the depth of coverage (DoC), heterozygote balance (Hb) and sequence coverage ratios (SCRs), as well as different forensic parameters of STRs and SNPs in a population study, were calculated to evaluate the overall performance of this new panel and its practicality in forensic application. In general, except for two STRs (DYS505 and DYS449) and one SNP (rs4288409) that performed poorly, the other 85 STRs and 293 SNPs in our panel had good performance that could strengthen efficiency for human identification and paternity testing. In addition, discordant STR genotype results between those generated from capillary electrophoresis (CE) and from the MPS platform were clearly illustrated, and these results could be a useful reference for applying these particular non-CODIS STRs in forensic practice.

Large-scale sequencing and comparative analysis of oenological Saccharomyces cerevisiae strains supported by nanopore refinement of key genomes.

Saccharomyces cerevisiae has long been part of human activities related to the production of food and wine. The industrial demand for fermented beverages with well-defined and stable characteristics boosted the isolation and selection of strains conferring a distinctive aroma profile to the final product. To uncover variants characterizing oenological strains, the sequencing of 65 new S. cerevisiae isolates, and the comparison with other 503 publicly available genomes were performed. A hybrid approach based on short Illumina and long Oxford Nanopore reads allowed the in-depth investigation of eleven genomes and the identification of putative laterally transferred regions and structural variants. A comparative analysis between clusters of strains belonging to different datasets allowed the identification of novel relevant genetic features including single nucleotide polymorphisms, insertions and structural variants. Detection of oenological single nucleotide variants shed light on the existence of different levels of modulation for the mevalonate pathway relevant for the biosynthesis of aromatic compounds.

BC-store: A program for MGISEQ barcode sets analysis.

Here we present the devised BC-store-a program for analyzing and selecting sets of barcodes for sequencing on platforms manufactured by MGI Tech (China). The app is available as an open source in Python3 and as a desktop version. The application allows analyzing the compatibility of barcodes on a single lane of a flow cell in a set in the case of equal and arbitrary fractions. In addition, with the help of this tool barcodes can be added to an existing set with custom share options. In this paper we describe how BC-store works for different tasks and consider the effectiveness of using BC-store in sequence lab routine tasks.